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Abstract 



We identify the leading term describing the behavior at large dis- 
tances of the steady state solutions of the Navier-Stokes equations in 
^ I 3D exterior domains with vanishing velocity at the spatial infinity. 

o : 1 Introduction 

We consider the 3D steady-state Navier-Stokes equations in exterior domains 
Q I and study the behavior of the solutions "near infinity" . The equations are 

X' divn = f ^^^ \^^- ^^-^^ 

where Bji^ denotes the ball of radius Rq centered at the origin. Our main 
assumption about the solutions will be the decay condition 

l^^^)l^^-T^ mM-n^Ro (1-2) 

Kq + \x\ 

for sufficiently small C*. The specific boundary conditions at OBrq will play 
no role in our results. 
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Naively one might think that the behavior near infinity of the above 
solutions should be given by the linearized equation. An immediate well- 
known objectioro to that is that we expect decay \V''u\ = 0(|a;|~'^~"'^) and 
{V^pI = 0{\x\~''~^) as X — > oo, and therefore the non-linear term in the 
equation should have the same order of magnitude as the linear terms, mak- 
ing the accuracy of the linearization questionable. This heuristics is made 
rigorous in [S], where it is proved that the leading order term describing the 
behavior of the solutions cannot be given by the linearized equation. 

In this paper we identify explicitly the leading order behavior of the above 
solutions near infinity. We show that it is given by the explicit solutions cal- 
culated by L. D. Landau in 1943, see [Sill]- These calculations were revisited 
and certain extensions were obtained in [U [TB] . The Landau solutions were 
recently characterized in [13] as the only solutions of the steady Navier-Stokes 
equation in M^ \ {0} which are smooth and (— l)-homogeneous in M.^ \ {0}. 
The Landau solutions in M.^ \ {0} can be parametrized by vectors 6 G M^ in 
the following way: For each fe G M'^ there exists a unique (— l)-homogeneous 
solution U^ of the steady Navier-Stokes equations together with an associ- 
ated pressure P^ which is (— 2)-homogeneous, such that U^,P^ are smooth 
in M^ \ {0}, U^ is weakly div-free across the origin and satisfies 

- AU'' + div([/^ ® U^) + VF^ = bd{x) (1.3) 

in the sense of distributions. Here d{x) denotes the Dirac function and we 
use the standard notation u^v ioi the tensor field UiVj defined by the tensor 
product of the vector fields u = {ui, . . . , Un) and v = {vi, . . . , f„). We also 
use the standard notation divT for the vector field -^Tij. (The uniqueness 

of U'^, P^ is much easier to prove if we add the requirement that U^ be axi- 
symmetric with respect to the axis passing through the origin in the direction 
of the vector 6, but as was shown in [13] , this additional symmetry assump- 
tion is not necessary.) As noticed by Landau, the solutions can be calcu- 
lated explicitly in terms of elementary functions, see formulae are (I3.1|) . (l3.4j) 
and (13.51) . (It was observed in [13] that the Landau solutions are in a natural 
one-to-one correspondence with the group of conformal transformations of 
the two-dimensional sphere, and Landau's formulae can be also derived from 
this observation by using some standard geometry.) 



^The objection was probably raised already by the classics in the 19th century. 



If u,p is a solution of fll.ip . we will denote by 

Tij = T,j{u, p) = p6ij + UiUj - (^ + Q^) (1-4) 

the momentum flux density tensor in the fluid. Our main result is the fol- 
lowing: 

Theorem 1. For each a G (1,2) there exists e = e{a) > such that the 
following statement holds true: Let u,p be a solution of fll.ip in M^ \ Bptg 
satisfying condition (11. 2p with C* < e. Let b = (61,62,^3) be defined by 

Tij{u,p)nj{x) (1.5) 

for some Ri > Rq. (Note that the integral is independent of Ri.) Let U^ be 
the Landau solution corresponding to the vector b. Then 

u{x) = U\x) + 0{\x\''') as \x\-^ 00 (1.6) 

and, for a suitable constant po, 

p(x) - po = P\x) + 0{\x\-''-^) as \x\-^ 00. (1.7) 

Remark 1. Standard estimates for the linear Stokes system (such as esti- 
mate (12. 5p in the next section), together with the scaling symmetry u{x) — >■ 
Xu{Xx) of Navier-Stokes can be used to show that any solution of (II. ip sat- 
isfying (II. 6p will also satisfy 

V^m(x) = V'=[/''(x)+0(|x|-'^-") as |x| ^00, for A; = 1,2,..., (1.8) 
and 

VV = V'T'' + 0{\x\-''-^-'') as \x\ ^ 00, for A; = 1, 2, . . . . (1.9) 
See for example [H] for an argument of this type. 

The existence of expansions similar to (11.61) with U^ replaced by a less 
specific term, namely a (—1)— homogeneous function, was is studied in [11]. 
The main result of that paper is, roughly speaking, that under smallness 



conditions similar to fll.2p . the solutions of fll.ll) are "asymptotically (— 1)- 
homogeneous" . As is shown in ^Jj, the leading term of the asymptoti- 
cal expansion at cxo of any solution of (11.11) which is asymptotically (— 1)- 
homogeneous must be given by a Landau solution. (This result remains true 
even for large data, since the proof is not based on perturbative arguments.) 
Therefore the results in [TT] together with the results in [13] imply a version 
of Theorem [H Our proof in this paper is much simpler than the proof one 
could get by combining [TT] and [13] . Also, if one tried to evaluate explicitly 
the values of the constants in the smallness conditions, the constants coming 
from the proof here would probably be more favorable. 

The proof of the Theorem [T] is given in the following sections. The main 
idea of the proof is as follows. Assume for simplicity that u satisfies the "no 
outflow to infinity condition" 

u{x)-n{x)=0 (1.10) 

for some Ri > Rq, where n{x) denotes the unit normal to dBji^. (Since u is 
div-free, the last integral does not depend on Ri.) We extend the fields u,p 
to fields defined in all M^ and satisfying the inhomogeneous equation 

divM = J 

The extension needs to be done in a way which enables one to control smooth 
norms of the extended function by the corresponding norms of the original 
function. We have b = Jj^a /. We then search for solutions of (11.111) in the 
form u = U'^ + v, where U'' is a suitable regularization of the Landau solution 
[/*. The equation for v is solved (for small data) by a standard perturbation 
analysis in the space of continuous functions with decay 0(1x1"°^) as \x\ — >• cxo. 
For this argument one needs both an existence result (for v) and a uniqueness 
result (for u), and therefore it looks unlikely that our method could be used 
in a large data situation. 

The general situation when the outflow J^^ u ■ n does not vanish can 
be handled by a standard method of writing u = a + v where v has no 
outflow, and a is a suitable multiple of the canonical outflow field 7%-. See 
for example [6], Section 2.2, [7], Chapter IX, or [TT], Remark 3.2. Roughly 
speaking, the part of the flow which produces a non-zero outflow has decay 



0(|x|^^), and therefore it does not influence tlie main term in fll.6p at large 
distances. See Section H] for details. 

One can see easily by looking at the Landau solutions that the best pos- 
sible decay rate a for which the result might still be true is a = 2. We 
conjecture that the result indeed remains true for a = 2. However, as the 
example in Remark [2] shows, one would probably need to go beyond the 
elementary perturbation theory used in this paper to prove that. 

The problem of steady-state solutions of the Navier-Stokes equations in 
exterior 3D domains has a long history going back to Leray's paper [TO] . 
Leray proved the existence of solutions with finite energy. Such solutions 
are easily seen to be smooth since the steady state equation is subcritical 
with respect to the energy estimate. However, the precise behavior of these 
solutions as |a;| — > oo is a more subtle problem. This problem shares some 
features with the (super-critical) regularity problem for the time-dependent 
equation, since there seems to be some vague duality between regularity 
(or short-distance behavior) of super-critical problems and asjTiiptotics at 
large times/distances (or long distance behavior) of sub-critical problems. 
In particular, in both cases it seems to be important to obtain some local 
control of the energy flux which is stronger than what one can immediately 
get from the known conservation laws. In this paper we will not address these 
difficult issues, which arise for large data, and we will only treat the small 
data situation, which can be handled by a simple perturbation theory, and 
is independent of the energy methods. For the steady state exterior problem 
this approach was pioneered by Finn, see e. g. P, [7]. We note that the 3D 
exterior problem with non-zero velocity at cxo has been more or less fully 
solved, even for large data, see [2], [7], since the non-zero velocity at infinity 
sufficiently regularizes the flow. In the 2D situation many problems remain 
open even in the case of non-zero velocity at inflnity, see [H [7] . 

How reasonable are our assumptions? We note that Finn proved in [6] the 
existence of solutions satisfying our assumptions for quite general boundary- 
value problems in exterior domains under smallness assumptions on the data. 
See also [7] for extensions of these results (still under smallness assumptions). 
It is quite conceivable that Theorem [T] remains true even without assuming 
that the constant C* in (11.21) is small. The proof of such a result would how- 
ever require to go beyond the perturbation theory and the standard energy 
methods. 

The following example, taken from [12], shows that the question of re- 
laxing the decay condition (II. 2p to a slower decay might be quite subtle. 



Consider the equation 

a 1 

- Am+ (1 -a)wVM+ -V|itp + -MdivM = in M^ , (1.12) 

for vector fields u in M.^. The number a G (0, 1) is a parameter. (For a = | the 
non-linear term in (11.121) can be written as div Q{u, u) for a suitable quadratic 
expression Q.) The equation has the same energy estimate as the Navier- 
Stokes equations. It turns our that (I1.12p has a nontrivial global smooth 
solution u satisfying \u{x)\ ~ |x|~^/^ as x — >• oo. Since it appears that the 
various perturbation and energy methods used for the steady Navier-Stokes 
should also work for (11.121) . at least in the case a = |, the properties of u 
indicate some limitations to these methods. On the other hand, we should 
remark that steady Navier-Stokes does have some special properties which 
are probably not shared by (I1.12p . For example, for steady Navier-Stokes 
the quantity \\u\^ +p satisfies a maximum principle, see e. g. pj. 

The paper is organized as follows: In Section [2] we explain the material 
necessary for extending the solutions to M? in a controlled manner. In Sec- 
tion [3] we recall the necessary facts about Landau's solutions. Finally, in 
Section m we explain the perturbation argument. 

Acknowledgement. The authors are indebted to Professor G. P. Galdi for his 
valuable comments on a preliminary version of the manuscript, which lead 
to significant improvements. 

2 Preliminaries 

We consider the solutions of the steady Navier-Stokes equation 

- Ati + uVu + Vp = , , . 

divM = ^ ^ ^ 

which are defined "in the neighborhood of infinity", i. e. in the region M^\S^g, 
where B^^^ denotes the ball of radius Rq centered at the origin. In this section 
we will be interested in the solutions which satisfy 

l^^^)l^^-T^ inM3\^ (2.2) 

Rq + \x\ 



and the "no outflow to infinity" condition 

u-n = (2.3) 



/ 



for some Ri > Rq. (We note that the integral in (I2.3P is independent of Ri, 
since divu = 0.) 

The constant C* above will play a special role and we distinguish it from 
the "generic constants" which will be denoted by c. If c depends on a pa- 
rameter X and we want to emphasize this dependence, we will write c{X) 
instead of c. The value of c can change from line to line. 

By a solution of (12. ip we mean a smooth function vector field u in M.^\Bji^ 
which satisfies (12. ip for a suitable p. (The pressure p will be considered only 
as a "secondary" variable: Instead of saying "the solution {u,p)" , we can 
just say "the solution u" , with the understanding that (12.10 is satisfied for a 
suitable p.) Various other notions of solutions are used in the literature (e. 
g. weak solutions), but under the assumption (12.20 they all coincide are are 
equivalent to the one defined above. 

One reason that the above way of thinking of p only as an auxiliary 
variable works quite well is that the linear steady-state Stokes system 



—Au + Vp = div / 
div -u = 

satisfies local elliptic estimates of the form 



(2.4) 



l|VM|U(B.o,fl) + lb- (P)i?.0.«ll^(-B.0,«) - 

c(i?,X)||/|U(ij^^,,,) +g(i?,X)||n||ii(5^^_,,) , (2.5) 

where B^q^r denotes the ball of radius R centered at xq, {j))b^ r is the 
average of p over the ball B^g^R and X can be any space in which classical 
elliptic estimates work, such as an L^-space with p G (1, oo) or a Holder space. 
The main point of estimate (12. 5p is that there is no p on the right-hand side. 
See for example [2] for details. 

The linear estimate (12. 5p combined with the standard bootstrapping and 
scaling arguments (using the scaling symmetry u{x) — > Xu{Xx)) imply that 
solutions of (12. ip satisfying estimate (12.20 with C:^ < M also satisfy 

\V''u{x)\<c{k,M) f*.,+, mM.^\B2R,ioTk = l,2,... . (2.6) 



We now relate the solutions of fl2.ll) in R^ \ Bfi^ to the solutions of the 
equation in M^ with non-trivial right-hand side: 



-Au + uVu + Vjo = / , 
div -u = 



m 



R\ (2.7) 



Let w be a solution of 1^^ in M^ \ Br^ satisfying (Q with C^ < M and 
let p be the associated pressure, defined up to a constant. Using (12. 6 p we see 
that we can in fact choose a "normalized" p so that, for C^, < M, we have 

I VVl < c{k, M)-—^^—— in M3 \ 5 for A; = 0, 1, 2, ... . (2.8) 
[Ro + 1x1)*-+"^ 

We can now extend u,p from M^ \ B^r^ to u,p defined in M^ such that 
div -u = in M^ and 

\V''uix)\<c{k,M)——^^-—- inM3forA; = 0,l,2,... , (2.9) 

together with 

\V'p{x)\<c{k,M) f*,,,, in R3 for A; = 0,1, 2,... . (2.10) 

The construction of the extension p — »• p is standard. To be able to con- 
struct the extension w — > it, we of course need condition (12. 3p . With (12.31) 
satisfied, the existence of a smooth div-free extension u (not necessarily sat- 
isfying (12. 9p ) is also classical. The construction of a div-free extension satis- 
fying (12. 9p can be carried out in many ways. One can proceed for example as 
follows: Let rj: [0,cx3) -^ [0, 1] be a smooth function such that r]{r) = for 
r < 2 and rji^r) = 1 for r > |, and let ijRoir) = vi^)- Now set it = rjR^^u + v, 
where w is a suitable solution of the equation div v = —uVriJi^ = g which is 
compactly supported in B^r^. The equation divf = g has of course many 
compactly supported solutions, but it is possible to construct a solution op- 
erator S: g ^ V = Sg which has the required regularity properties. (Such 
an operator is sometimes called a Bogovskii operator.) See for example [3] 
or [7], Chapter III. 3 for details. 

We have 

-Au + uVu + Vp = / , 

div -u = 
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I inM^^ (2.11) 



where the right-hand side / is supported in B^ji^ and satisfies 

|V^/(x)|<c(A:,M) ^* in M3 for A; = 0,1, 2,... . (2.12) 

(i?0 + IXl)"-^-^ 

Dropping the tildes and changing Rq, if necessary, we see that the study of 
solutions of Navier-Stokes defined in the neighborhood of infinity and satis- 
fying the "no outfiow" condition (12.31) and the growth condition (12.21) can be 
reduced to the study of the solutions u of the inhomogeneous equation (12.71) 
with / supported in Br^ and satisfying (12.121) . and u satisfying (12.21) globally, 
with C* replaced by cC^. 



3 The Landau Solutions 

The Landau solutions are smooth (— l)-homogeneous solutions of the steady- 
state Navier-Stokes equations defined in M'^ \ {0}. Under the additional 
assumption of axial symmetry, these were first calculated by L.D. Landau in 
1943, see [H [9]. In |I3] it was proved that we do not get any new solutions 
if the assumption of axial symmetry is dropped. To write down the explicit 
formulae, we will use the standard polar coordinates r, 9, (p defined by 

Xi = r sin 9 cos (f , 
X2 = r sin 9 sin ip , 
X3 = r cos 9 . 

The explicit formulae in polar coordinates for the Landau solution U and 
the corresponding pressure P are as follows: 

2 sine 



Ue 
P 



r(A-cos9) ' (311 

t/^ = , ^ • ^ 

4(Acos6'-l) 
r2(A-cos0)2 ■ 



In the above formulae, A is a parameter satisfying A> 1. The velocity field 
U can also be expressed in terms of the stream function 

2r sin^ 9 



as 

IJ = 1 9ip 

^ r sin 8 rdd ' 

Ue = m^ , (3.3) 

" rsmO or ' ^ ' 

U^ = . 

The integral curves of the velocity field U are given the the equations tp = 
const, and ip = const. 

Clearly U, U ®U and P are locally integrable, and a direct calculation 
(see e. g. j9]) gives 

- At/ + div(f/ o [/) + VF = p{A)e:i5{x) , (3.4) 

where 63 is the unit vector in the positive xs-direction and 

/ 1 4 — 1 4/1 \ 

^(A,^16.(^ + 5^^.og— + 5^^^^). (3.5) 

It is not hard to check that the function (3 {A) is monotonically decreasing 
in (1, cxd) and maps this interval onto (0, 00). In particular, (3 has an inverse 
function 7: (0, 00) ^ (1, cxo). 

It is instructive to compare the formula (13. 2p with the corresponding 
formula for the linear Stokes system. Namely, the solution of 

-At/nn + VPHn = e35(x),| .^^3 (g.g) 

divf/iin = J 

satisfying ?7iin(a^) — > as x — > 00 is given by the stream function 

^lin = -— rsin^6' . (3.7) 

Stt 

We can get another useful comparison if we express the solution of the prob- 
lem 



,. „ ^ inM3, (3.8) 

with the condition Us{x) — > as x ^ 00 in terms of the Landau solutions. 
The formula (for e > 0) is 

1 , , , 

Ue = -U\A='y{e) (3.9) 

10 



and -Ug is given by the stream function 

r sin^ 6 

^e = —T-. -, ■ 3.10 

£7(£) — 5 cost/ 

We will now regularize the Landau solutions near the origin in the fol- 
lowing way. Let tq > 0. Consider a smooth function p: [0, oo) — > [0, cx)) such 
that p{r) = for r < Tq, p{r) = r for r > 2ro and the A;— th derivative p'^^\r) 
is bounded by c{k)r^~'^, and define 

~ ~ 2p(r)sin^g 

^ = ^-^'-°= A-cos^ • ^'-''^ 

With the help of -0 we now define the regularized velocity field f/ = UA,ro by 
the formulae (13. 3p . with -0 replaced by -0. We also define 

p_p _ 4p(r)(Acosg-l) 

^"^^•'■"" r3(A-cos^)2 • ^^-^^^ 

It is easy to check that for A > Aq > 1 we have 

\^'U\< ./^^'mL mM3forfc = 0,l,... (3.13) 

and 

l^'^l ^ ./^^'mL i^ M^ for fc = 0, 1, . . . . (3.14) 

So far we have mostly considered the Landau solution which are axi- 
symmetric with respect to the x^ axis. However, it is clear from the above 
that for each non-zero vector 6 G M^ there exist a unique Landau solution U^ 
and the associated pressure P^ which are axi-symmetric with respect to the 
axis M • h and satisfy 

-AU^ + div([/^ ® U^) + VP^ = b 6{x) , 

div U^ = . 



b _ n (3-15) 



We also set U^ = 0. For each [/*, P* the above construction of the regularized 

/,^„ and P' = P,'^ 



solutions gives the regularized fields U^ = U^ and P^ = P^l which, for 



|6| < M will satisfy the estimates 

^ 

[ro + \x\Y+^ 



V^U^\<c{k,M)- ^^T-T inM^for A; = 0,1,... (3.16) 



and 



jkBbi / „/,. .^^ 1^1 



IV'^P"! <c(A:,M)- ',' ^^, inM3forfc = 0,l,... . (3.17) 

(ro + |x|)''+^ 



11 



4 Perturbation Analysis 

Let / be a sufficiently regular compactly supported vector field in M^. Let 
b = /ig3 /. Let ro = 1 and let U = U^ = U^^ and P = P^ = P^^ be the 
regularizations of the Landau solutions U'', P^ corresponding to the vector 
b constructed in the previous section. We will seek solutions of the steady 
Navier-Stokes equation 

-Au + vS/u + Vp = / 



n f inM^ (4.1) 

in the form u = U + v. We set 

F = -AU + UVU + VP . (4.2) 

It is easy to check that Jjgg F = b = f^^ f. The equation for v becomes 

-Av + UVv + vVU + vVv + Vq = f-F 



f in M^ . (4.3) 

divf = ' 

Let us choose a fixed a G (1,2). We will prove that under some smallness 
assumptions equation (14.31) has a unique solution v with decay 0(|x|^") as 
|a:;| —>■ oo. An important point is that, by our construction, Jj^aif — F) = 0. 
Using the scaling symmetry, we see that we can assume tq = 1 without loss 
of generality. 

Let G = Gij be the Green tensor of the linear Stokes operator. We 
note that the vector field Gjs is given by the stream function (13.71) . Another 
explicit formula for G is 

G,W = i(-.,A + ^)N (4.4) 

For our purposes here we will only need the following obvious estimate 

|VG(x)| < ^ . (4.5) 

The required solutions of (14. 3 p will be found for small data by a standard 
perturbation argument. Let X^ be the space of all continuous div-free vector 
fields u in M^ satisfying u{x) = 0(|a;|~") as x — > oo. A natural norm in X^ 
is given for example by 

\u\a = snp {1 + \x\)'' \u{x)\ . (4.6) 

X 

Our perturbation analysis is based on the following elementary estimates: 

12 



Lemma 1. Using the notation above, let b = J^^ f be the vector used in the 
construction ofU. Then for \b\ < M we have 

[G *diY {U^v + v^U)U<c{a,M)\b\[v\^ , (4.7) 

and 

[G * div(t; ® w)U < c{a)lvUwU • (4.8) 

Proof. The proof these estimates is standard. We move the derivatives to 
G, use the definition of the norm, after which the only remaining task is to 
estimate the integral 

dy 

\x-y\'^{l + \y\ 



nx) = I ,^ ..„./, ,..,..■,, , (4.9) 



where /3 G {1, a}. It is enough to consider only the case P = 1. Clearly I{x) 
is bounded for \x\ < 1. To estimate I{x) when |x| is large, let us write x = te 
with |e| = 1 and make the substitution y = tz in (14.91) (with P = 1). We 
obtain 

'('») - '-" L ie-.p(/..|3ir. ^ '- L ^3^ ■ (^-10) 

Since we assume a G (1,2), the last integral is bounded, and we see that 

t"J(te) < c{a) . (4.11) 

Combining this estimate with the estimate of I{x) for |x| < 1 we see that 

(1 + |x|)"/(a;) < c(a) for all x eR^ . (4.12) 

This completes the proof of estimates (14.71) and (14. Sp . D 

We have shown that the linear operator 

Tq: V ^G* div([/ ®v + v^U) 

is continuous from Xa to X^, and its norm is bounded by c(a, M)|6|. Also, 
we have shown that the bi-linear operator 

B: {v,w) ^ G * div(t; (g) w) 
13 



is continuous from X^ x X^ — > Xa, with the bound 

lB{v,w)\^<c{a)lvUwU ■ 
We let y = G * (/ — F) and re-write equation fl4.3p as 

v + Tq{v) + B{v,v) = V . (4.13) 

Since J^sif — F) = 0, we have V = 0(|2;|~^) as |x| -^ oo. Standard pertur- 
bation arguments (such as the Imphcit Function Theorem) now imply that 
equation (14.131) has a solution v when V is sufficiently small in Xa- (A sim- 
ple sufficient condition for that is that, in addition to J^sif — F) = and 
the restriction on the support on / — F, the field f — F he small in L2~^^ 
with some 6 > 0.) Moreover, the solution is unique in some small ball in Xa 
(centered at the origin). These statements can be made more quantitative 
if we use the special form of the perturbation (namely that it is quadratic 
in v). For example, one can use the following folklore lemma: 

Lemma 2. Let X be a Banach space. Let T: X ^ X be linear with \\Tx\\ < 
e\\x\\ for all X G X, and let B: X x X ^ X be bilinear with \\B{xi,X2)\\ < 
c|fi|| ||x2|| for all a;i,X2 G X. Let y ^ X with \\y\\ < ^J . Let < 
G < C2 be the two roots of the equation ^ = \\y\\ + £^ + ci^ , i- e. ^1^2 = 

^^^—2^ ■ Then the equation 

x + Tx + B{x,x) = y (4.14) 

has a solution x satisfying \\x\\ < ^1. Moreover, the solution x is unique in 
the open ball {x E X,\\x\\ < ^2}- 

Proof. The proof is standard and we include it for the convenience of the 
reader. Consider the map F(x) = y — Tx — B{x,x). We have ||F(x)|| < ||?/|| + 
£||x|| + c||a;|p which shows that for ^1 < ||x|| < ^2 we have ||F(x)|| < ||2;|| 
and that, for any 6 > 0, the iterates F{x),F'^{x) = F{F{x)), . . . , F^{x), . . . 
enter the ball of radius ^i + d after finitely many steps. At the same time, 
we have ||F(xi) — F(x2)|| < £\\xi — X2II + c\\xi — a;2||(||a;i|| + ||x2||) which 
shows that F is a contraction of any closed ball of radius ^ G [^1, 2 )- '-' 

Proof of Theorem 1. Let us ffist assume that the "no outflow to infinity" 
condition (11.101) is satisfied. In this case the statement of Theorem [1] is a 
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direct consequence of the construction of the extensions in Section [2] and 
Lemmata [T] and [2J Note that we not only need existence and uniqueness 
for V in (14. 3p . but we also need uniqueness (with smallness assumptions) 
for u in (14. ip . The uniqueness of u in our situation is well-known (see e. 
g. [6] or [7J), and can also be easily proved from Lemma [2] and (an obvious 
modification of) Lemma [H 

The situation when we have some outflow to infinity can be handled by a 
standard method of using the canonical outflow field ^, see for example [6], 
Section 2.2, [7], Chapter IX, or [llj. Remark 3.2. Assume Rq = 1 without 
loss of generality. Let a be the multiple of the vector field -A^ which has 
the same outflow as u. Note that a satisfies the Navier-Stokes equation (II. ip 
in M^ \ {0} with the associated pressure field tTq = — ||ap. Let us write 
u = a + w and p = Tfa + Pw The field w satisfies the no outflow condition, 
and we can extend it to a div-free field w with the control similar to (12. 9p . 
We can also regularize a and tTq in Bi (while not changing them outside Bi) 
so that estimates similar to (12. 9p and (I2.10p are satisfied. Let us denote be a 
and tTq these regularized functions. Finally, we extend pw to pw with control 
similar to (I2.10p . Let u = a + w and p = T^a + Pw (Note that u is not 
div-free in Bi.) Let / = divT, where T = T{u,p) is given by (ll.4p with 
u,p replaced by u,p. We note that the vector b given by (II. 5p can also be 
expressed as 6 = f^^ f . We will now search a div-free vector field z and a 
function p^ satisfying div T(a + z,7ra + Pz) = /• We seek z,pz in the form 
z = U'' + v, Pz = P^ + q, where [/'', P^ are the regularizations of the Landau 
solutions U^, P^ constructed in Section [31 It is now easy to check that the 
perturbation theory of Section H] gives the required solution. 

D 

Remark 2. The borderline space Xa in which a more sophisticated per- 
turbation analysis might possibly work is the space X2. (This corresponds 
to the naturally expected decay 0(|a;|^^) for v.) However, a perturbation 
analysis in X2 cannot be based only on the decay properties of U (as was 
the case with our simpler analysis for a < 2). To see this, let e G (0, 1) and 
consider the equation 

si Tl — S) X 

-Au+- -dw{-—^u) = f{x) inR". (4.15) 

One can check by direct calculation that when / = the function xilxl"""*"^ 
is a solution of this equation away from the the origin. Let 77 be a smooth 
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function in M" which vanishes in the unit ball and is equal to 1 outside of the 
ball of radius 2. An easy calculation shows that the function u = rj Xilxl'"-^^ 
satisfies fl4.15p with J^^ / = . Moreover, one can change the coefficients of 
the equation in the unit ball so that they become smooth. 

To get results in the space X2, one would probably have to prove optimal 
decay estimates for the linear equation 

-A. + ^V. + vVU + Vg = f-F,) .^ ^3 (4.15) 

divv = J 

by a non-perturbative approach, and then treat the quadratic term in (14.31) 
perturbatively. The above example shows that to get the optimal decay 
0(|x|~^), one would need to use more information about U than just its 
decay properties at cxd. We conjecture that for large |x| the perturbation v 
from the Landau solution U^ indeed has the decay v{x) = 0(|x|^^), at least 
for small data. 
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